epl draft 



Fractal space frames and metamaterials for high mechanical effi- 
ciency 



o 
in 

(N 



o 

CZ2 



R. S. Farr 1 and Y. Mao 2 

Unilever R&D, Olivier van Noortlaan 120, AT3133, Vlaardingen, The Netherlands 
School of Physics and Astronomy, University of Nottingham, Nottingham, United Kingdom 



-a 
c 
o 

(N 
> 
O 

o 
o 



PACS 46 . 32 . +x - Static buckling and instability 
PACS 46 . 25 . Cc - Static elasticity: theoretical studies 
PACS 46.70.Hg - Membranes, rods and strings 

Abstract. - A solid slender beam of length L, made from a material of Young's modulus Y 
and subject to a gentle compressive force F, requires a volume of material proportional to L 3 f 1 ^ 2 
[where / = F/(YL 2 ) <C 1] in order to be stable against Euler buckling. By constructing a 
hierarchical space frame, we are able to systematically change the scaling of required material 
with / so that it is proportional to L 3 y T ( G + 1 )/( G + 2 ) ( through changing the number of hierarchical 
levels G present in the structure. Based on simple choices for the geometry of the space frames, 
we provide expressions specifying in detail the optimal structures (in this class) for different values 
of the loading parameter /. These structures may then be used to create effective materials which 
are elastically isotropic and have the combination of low density and high crush strength. Such a 
material could be used to make light-weight components of arbitrary shape. 
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j | Introduction. — It has recently been shown that frac- 
tal [TJ design principles can be used in structures to im- 
prove their mechanical efficiency under conditions of gen- 
tle compressive loading [2j[3]- The purpose of this paper 
is to introduce much simpler structures which display the 
same phenomenon, are more practical to manufacture, and 
which can then be used as components to construct an 
effective composite material (metamaterial) . This mate- 
rial will, on long length scales, be elastically isotropic, and 
have a very low density but a comparatively high compres- 
sive strength or 'crush pressure' 0]. It could in principle 
then be used to manufacture light-weight components of 
any 3d form, in the same way that solid foams are used cur- 
rently [5]. Being composed of connected, straight beams, 
this structure should be amenable to modern fabrication 
techniques [BHE]- 

To begin the argument, we consider fractal trusses or 
space- frames 0] . We wish to design a structure made from 
a volume V of material of Young's modulus Y, which is 
required to support a compressive force F applied at freely 
hinged end-points, separated by a distance L. It is then 
useful to define a non-dimensionalized force parameter 



and non-dimensionalized volume 



V/L 3 < 1. 



(2) 



The aim is to find a structure which, for small values of /, 
minimizes v. In particular, the scaling of achievable values 
of v with /, as / —> is of interest. 

For a solid cylindrical column, the limitation is Euler 
buckling [9], which leads to v = 2n~ 1 / 2 f 1 / 2 (in contrast to 
the much higher efficiencies of v oc / possible for structures 
under tension) [HHU]. 

A space-frame. — Instead of a single solid cylindrical 
column, let us consider a space-frame. For ease of analy- 
sis (rather than because this is the global optimum) the 
structure is composed of beams of equal unstressed length 
Lq, which are freely hinged to each other at their ends. 

In cartesian co-ordinates, the end points of the structure 
lie at the origin and the point (0,0, L). The geometry of 
the space frame consists conceptually of a regular tetra- 
hedron at each of its ends, and between these a stack of 
n > regular octahedra, as shown stereographically in 
fig- DI a ) an d (b) for the cases n = and 3 respectively. 
Therefore the length of each component beam is 



/ = F/(YL 2 ) « 1 



(1) 



(3) 
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Fig. 1: Stereograms of 'generation F space-frames, with (a) 
=0 and (b) ni t i —3. To view these figures, hold the page 
about 30cm away, and look through the page, until the two 
images merge. 

A consequence of the rigid polyhedron theorem [TT] is 
that space frames which are convex triangular polyhedra 
composed of rigid beams are themselves rigid, and there- 
fore this design based on tetrahedra and octahedra ensures 
that there are no soft modes for the structure: any flexibil- 
ity must involve stretching or compressing the constituent 
beams. 

When a small compressive load F is applied to the ends 
of the space frame, some of the component beams will be 
under tension, and others under compression. In particu- 
lar, the beams which are parallel to the x — y plane are all 
under tension, and the rest are under compression. The 
six beams connected directly to the ends are all under a 
compressional force 



-Fcom = F/V6, 



(4) 



while the other beams under compression support only 
half this load apiece. 

For the beams under tension, if n = 0, then all three 
are subject to a force of 2F/(3\/6)- On the other hand, 
if 72 > 1, the six tension beams closest to the end points 
are subject to a force F/(2\/6), and if n > 2, then all the 
remaining tension beams are subject to a force F/(3\/6). 

Each of the beams will behave as a linear spring under 
stretching or compressional forces small enough to avoid 
Euler buckling of these individual beams. We now sup- 
pose that all of the beams have the same effective spring 
constant (force per unit extension), given by k. Once more 
this choice is made for simplicity of analysis: we note as 
an aside that more detailed calculation shows that a gain 
in efficiency by a numerical factor of order unity is possi- 
ble by choosing some beams to have different spring con- 




Fig. 2: Space frame consisting of 36 octahedra composed of 
equal beams, connected to form a circle, and allowed to relax. 
The elastic energies for space frames of this type (and vary- 
ing numbers of octahedra) are used to estimate the bending 
stiffness of the space frame, treated as an equivalent slender 
beam. 



stants; however, the functional dependence of v of / is not 
affected. 

With this assumption of equal spring constants fc, cal- 
culation shows that the spring constant of the entire space 
frame under compression is given by 



K 



36fc 



lln + 43 



(5) 



provided n > 3. 

Next, provided n is large, the entire space frame will 
resemble a long, slender beam, with a bending stiffness 
YI, where / is the second moment of the area about the 
neutral axis for the equivalent slender beam [T^] . Because 
of the three-fold symmetry of the space frame on rota- 
tion about the z-axis, and the possible symmetries of the 
second moment calculation, we would expect the bending 
stiffness to not depend on the direction of bending. 

Finally, we would expect that as n — > oo, 



YI 



BLlk, 



(6) 



for some numerical constant B. This constant is diffi- 
cult to calculate analytically. However, by simulating long 
space frames which have their ends joined to form a circle, 
and which arc allowed to relax (fig. [2]), the relaxed energy 
for different numbers of constituent octahedra can be cal- 
culated. By choosing the number of octahedra up to 128 
and extrapolating to higher values, we find 



B = 0.245 ±0.001. 



(7) 



During Euler buckling of a freely hinged beam, the cur- 
vature at the end points will vanish |12j (this follows di- 
rectly from the freely hinged condition). Therefore using 



p-2 



Fractal space frames and metamatcrials 



eq. © in the expression for the maximum force F\, uc which 
a slender beam can sustain before buckling (T2"] : 



Fh 



n 2 YI 



(8) 



will provide a good approximation to the failure criterion 
for the whole space frame under a global buckling insta- 
bility, even when n is not large. 

Each component beam is also vulnerable to buckling 
(in fact the component beams which are joined to the 
end points are most vulnerable) , and therefore the force F 
which the structure can withstand is subject to two con- 
straints: the individual beams must not buckle, and the 
space frame itself should not undergo a global buckling in- 
stability. In order to be optimally efficient, the structure 
should be on the verge of both instabilities. 

Once we have chosen L, the space frames we have just 
described are specified by two parameters: the number of 
octahedra n, and the radius r of the circular cross-section 
of the component beams. This last is related to the spring 
constant of the beams through 



kL 
irY 



1/2 



(9) 



[with Lq coming from eq. ©], 

To design an optimal structure, one should therefore 
specify L, Y and F (or, more elegantly, the parameter 
/), and then choose the values of n and r which minimize 
v. This approach (especially for the hierarchical struc- 
tures described below where there are different n's), leads 
to equations which are difficult to solve analytically. In- 
stead, we construct optimal structures through a simpler 
algorithm, which is possible because of the minimal cou- 
pling between hierarchical levels in the structure. This 
generalizes easily to more complex structures, and is im- 
plemented as follows: 

First, we parameterize the optimum space frames 
through a new dimensionless parameter 



/o = (F/V6)/(YL 2 ), 



(10) 



which is the / parameter calculated for one of the single 
constituent beams which is most vulnerable to buckling. 

The condition that this beam be at the Euler buckling 
limit then allows us to calculate r, giving 



r = Lq 



4jo 



1/4 



(11) 



Second, the condition that the entire structure not 
buckle globally allows us to determine n, by using eqs. (|3"j). 
©and © - CtT]): 



3^ 
2 



B 1/2 f Q 



1/4 




-2+Ll.29/- 1/4 j. 



(12) 



Fig. 3: Generation 2 space frame, with 712,1 = 4 and 712,2 = 2. 

where ['J is the floor function. 

The third step is to calculate / and v in terms of /o, 
which follow directly from eqs. ©, © and (|T2"j) as 

/« 3.67/o (Ll.29/- 1/4 j)~ 2 , (13) 

^18.7(-l+Ll.29/- 1/4 j) (Ll.29/- 1/4 j)~ 5 . (14) 

For small /, we therefore have 

v^f 2/3 , (15) 

which represents a qualitative increase of efficiency over a 
solid cylindrical column, for which v oc f 1 ^ 2 . 

Hierarchical space frames. — In the last section, we 
described a simple space frame, consisting of solid cylin- 
drical beams. We now re- name this a 'generation G = V 
structure, which for fixed L can be specified by the ra- 
dius r of the cross section of the constituent beams, and 
the number of constituent octahedra, which we now re- 
name Tii,!.. In this notation, the first index refers to the 
generation number G, and the second to the depth in the 
hierarchical structure (see below). 

To define a generation G = 2 structure, we replace all 
the comprcssional beams in a generation 1 structure, by 
(scaled) generation 1 space frames, which are all identi- 
cal to each other. The number of octahedra present in 
the two hierarchical levels of the structure may however 
be different; with 712.1 being the number of octahedra in 
the smallest component space frames, and 712,2 being the 
number of large octahedra (composed in part of composite 
beams) in the entire structure. However, for simplicity of 
manufacture, the tension members are allowed to remain 
as simple solid cylinders, since no gain in efficiency is de- 
rived from making them more complex. An example is 
shown in fig. [31 with 712,1 = 4 and 77-2,2 = 2. 

This process may be repeated, in each iteration replac- 
ing the smallest comprcssional beams by identical (scaled) 
generation 1 space frames, and so obtaining structures 
with higher and higher generation number G. fig. [4]shows 
a generation 3 space frame, with 773^ = 4, 773 2 = 3 and 
^3,3 = 2. 

In the final generation G structure, all the comprcssional 
beams are solid cylinders of equal length Lq,o, and have 
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Fig. 4: Generation 3 space frame, with ri3,i = 4, 713,2 = 3 and 
^3,3 = 2. 



equal radii r. These are the smallest beams in the struc- 
ture. They have a bending stiffness YIq,o = 7rY> 4 /4, a 
spring constant before buckling of ko,o = nYr 2 / Lg,o and 
they support a maximum comprcssional force of Fg,o- 

The smallest beams make up small space frames which 
have a geometry similar to generation 1 structures. Each 
of these space frames are of length Lg,i, and contain nc,i 
octahedra. They have a bending stiffness YIg,i, spring 
constant kc,i and support a maximum compressive load 
of Fq \. 

Higher levels in the hierarchical structure are defined in 
the same way, so that eventually 



L = Lg,g 
F = Fg,g- 



(16) 
(17) 



From eqs. ©, (JB]) and ([5]), we can write down recurrence 
relations for all these quantities, valid for 1 < i < G: 



Lga = 



1/2 



(n G ,i + 2)L G ,i 



YI Gli = BL% i^kci-i 



36fc G ,»-i 
lln GiI + 43 

Fga = VsFga-i- 



(18) 
(19) 
(20) 

(21) 



There will also be tensional beams in the structure of 
various lengths, which have been 'left over' - i.e. not re- 
placed by composite beams - in the process of constructing 
the structure. The radii £g,; of these tensional beams are 
chosen (for simplicity of analysis) to give them the same 
spring constant as the comprcssional (usually composite) 
beams at the same level in the hierarchical structure. Thus 



the smallest tensional beams are identical to the compres- 
sional beams 

t G ,o = r, (22) 

and in general the tensional beams at level i £ (0, G — 1) 
have length Lga an d radius 



tGA 



irY 



1/2 



(23) 



Note: there is no tensional beam at level G, since this is 
the entire structure, which is under compression. 

Optimization. — As with the simple space frames 
above, we parameterize optimized space frames through 



fo = Fga/ (YLqq) 



(24) 



The first stage is to find the optimal choices for nc,i at 
each level in the structure. To do this, let us take Lg,o 
as fixed (we will determine it in terms of L only at the 
end of the calculation). We can then determine the radius 
of the smallest compression beams though imposing the 
condition that they be on the verge of Euler buckling: 



4/0 



1/4 



so that 



fc G ,o = 2Tr- 1 / 2 L G ,oYfo 



1/2 



(25) 



(26) 



The condition for buckling to not occur at level i £ (1, G) 
in the structure is 



Ft 



G.i 



< 



n 2 YI Gii 

T 2 ' 

GA 



(27) 



which from eq. (0, (JTg]) and (gUJ) leads for i = 1 to 



n-G,i 



-2+Ll.29/- 1/4 j, 



(28) 



and for 2 < i < G, we find from eqs. (|18H21|) . ([26)) and 
(j27|) that 



n G A = -2 



where 



An G ,j + 43 



a/o (2V3)* n (? 



£1/2^3/4 
A = —TTT7T- ~ 0.373. 

2(6 V4 



(29) 
(30) 



eqs. (|28p and (|2"9"|) allow us to calculate all the values for 
numbers of octahedra at different levels in the structure. 
A simple calculation (including the volume of the tension 
beams) then leads, for G > 2 to 



/07\ G / 2 G 

f=(i) ^ +2)- 



(31) 
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Table 1: Example calculation for the mass M if a structure 
required to support F = lOkN over a distance of L — 200m 
when the structure is made from a material similar to steel, 
with Y = 210GPa and density p = 8000kgm~ :! . This corre- 
sponds to / = 1.2 x 10~ 12 . 
G f n G ,i n Ga n G ,3 M 



TT 

6.63 x 1CT 9 
4.92 x 10~ 7 
6.36 x 1CT 6 



1.2 x 10 
1 

2 
3 



140 
46 
23 



47 
23 



24 



79 tonnes 
2920kg 
1790kg 
2180kg 



1 /243 

G-l ( q 

3+e n 

9=1 [j=l 



G/2 



n 



\ {n G ,j + 2) 3 



fo 



1/2 



2 2 «(n Gj + 2) 2 



(lln G>i + 43)(n GJ + l) 



(32) 



To illustrate these calculations, consider a space frame 
of length L = 200m which is required to support a force 
of F = lOkN, and which is made from a model mate- 
rial, similar to steel, with Y = 210GPa and a density of 
8000kgm~ 3 , so that / = 1.2 x 10~ 12 . 

A cable supporting this force under tension would re- 
quire a mass of 8kg (assuming a yield stress for the ma- 
terial of 200MPa, and neglecting the mass of couplings at 
the ends). The masses (M) of "steel" required for various 
structures described in this paper are shown in table [T] 

More generally, eqs. (pTTj) and (pi2"j) show that for fixed 
G, then as / — > 0, 



V oc / 



( G+ l)/(G+2) 



(33) 



Furthermore, for any particular /, there will be an op- 
timum generation number G op t, which gives the most ef- 
ficient (smallest v) structures in this class, table [2] shows 
the approximate ranges of / which correspond to different 
values of G op t- 

Finally, we note that there is still considerable opportu- 
nity to optimize the space frames presented here; the de- 
sign in this paper has been motivated by a search for the 
scaling behaviour in the limit of small /. If one wished to 
reduce the pre-factor, there are opportunities through al- 
lowing different beams to have different thicknesses in the 
structure (currently the structure has weak points near 
the ends), through choosing a different basic space frame 
geometry, and potentially also through moving from freely 
hinged to rigid boundary conditions at some of the joints. 
Such improvements will also also alter the crossover values 
in table H 

An elastic metamaterial. — For application to an 
elastically isotropic composite material, we note that for 
example, a face centred cubic lattice built from freely 
hinged beams between nearest and second neighbours [T3] 
has this property of isotropy. If we replace all these beams 
with the hierarchical structures described in this paper, 



Table 2: The approximate ranges of / and fo, for which dif- 
ferent generation numbers G op t give the global optimum space 
frames of the type described in this paper. 
Gopt range of / range of fo 



10 



</ 

10" 9 < / < 10~ 4 
10" 14 < f < 10" 9 
10" 62 < / < 10" 14 



io- 
io- 
io- 

10" 



</o 

< fo < 10~ 3 

< fo < 10~ E 
5 < fo < 10- 



then the resulting effective material will use a volume frac- 
tion of real material, where <f> is of order v. The crush 
pressure of the composite material will, from eq. (|33[) . scale 
with volume fraction according to 



Pc 



Ycj) 



(G+2)/(G+l) 



(34) 



For small enough <fi, the value of p c according to eq. ([341 
is much larger than a material built from a lattice of sim- 
ple beams. We note that cancellous bone (which is fractal 
[14] ) also has an unusual scaling of (anisotropic) compres- 
sive strength with density [5]- 

Other detailed considerations. We note that in 
the initial specification of the space frame, we have re- 
quired all the component beams to be freely hinged to 
one another. Potentially, these hinges require a certain 
amount of additional material. However, the extra mate- 
rial to make couplings for the ends of a beam is of order 
the cube of the beam diameter [TU]. Therefore, including 
all the small beams, we expect this consideration to in- 
crease the amount of material in the overall structure by 
a factor of (1+e), where e = ord(r/X Gj o) = ord(/d^ 4 ) <C 1, 
provided fo 1. 

The main reason for imposing freely hinged joints is that 
without this degree of freedom, the space frame would de- 
form under loads, leading to the beams being no longer 
straight, even before the nominal buckling load. Such de- 
formations could potentially weaken the structure. How- 
ever even with freely hinged joints the structure will de- 
form before failure, and although all the beams remain 
straight, the directions of force transfer will be slightly 
altered, which could disturb the calculations presented 
above. We would therefore like the elastic deformation 
of the structure at the time of failure to be small. Once 
again, we can check this directly: the compressional strain 
of the smallest beams at failure is [from eq. ([25])] propor- 

1/2 

tional to / <1, provided /o<l. 

Lastly, we make an observation on the geometry of the 
structure: if we have a space frame as described in the sec- 
tions above, with large G, and all the numbers n Gl , equal, 
then this structure will be a fractal over a suitable range of 
length scales. Furthermore, the Hausdorff dimension [15] 
will be a (decreasing) function of n G ,i- However, we see 
from eq. (|29p that for optimized structures, n G) i depends 
on the depth i in the hierarchy: in particular, as / — > 0, 
we find n Gj i oc (12/11)^ -1 ^ 2 . This leads to a change of 
Hausdorff dimension with length scale. 
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Conclusions. — We have presented an hierarchical 
design for space frames which allows us to systematically 
change the scaling of material needed versus compressive 
force under conditions of gentle loading. These can then 
be used to construct a light, strong mctamaterial. 

The resulting expressions for space frame efficiency are 
similar to those of curved shell struts [3], but the latter 
are more efficient than the designs presented in this pa- 
per. Considerable scope for further optimization of the 
pref actors however remains. 

Both these approaches are an attempt to solve an op- 
timization problem which is rather simply stated: 'given 
the loading condition specified by /, what geometry of 
an elastic material minimizes i>?'. On the basis of recent 
work in this area, the authors conjecture that for small 
enough /, the answer will consist of some kind of hier- 
archical structure. However, mathematical tools for ad- 
dressing optimization problems of this kind are notable 
mainly for their absence, and even recent numerical ap- 
proaches [T6lfT7] seem not to be powerful enough to tackle 
the problem directly. 

* * * 

The authors wish to thank Dr Joel Segal of Notting- 
ham University and Andrew Matthews of Devon County 
Council for useful discussions. 
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